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The reasons which restrict opportunities of classical mechanics at the description of nonequilibrium 
systems are discussed. The way of overcoming of the key restrictions is offered. This way is based 
on an opportunity of representation of nonequilibrium system as a set of equilibrium subsystems. 
The equation of motion and the general Lagrange, Hamilton and Liouville equations for subsystems 
have been obtained. The way of a substantiation of thermodynamics is offered. 
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I. INTRODUCTION 

The development of the theoretical bases of physics 
of open nonequilibrium systems is collided with contra- 
dictions between classical mechanics and thermodynam- 
ics. These contradictions are most brightly displayed in a 
problem of a substantiation of the second law of thermo- 
dynamics or, in other words, explanations of irreversibil- 
ity. Since L. Boltzmann and to this day the attempts to 
solve this problem are undertaken [1-4]. 

In the basis of the most widespread contemporary ex- 
planation of irreversibility the mixing property inherent 
at the Hamiltonian systems is used. The mixing leads 
to the irreversibility if one postulates averaging of phase 
space on physically small volume. But the explanation 
of the nature of such averaging within the framework of 
classical mechanics is impossible [2, 3]. 

With the purpose of search of solutions of these prob- 
lems the hard disks systems were investigated [5, 6]. Its 
studying revealed that if the nonequilibrium system of 
disks has been represented by interacting equilibrium 
subsystems (lES) then equilibration is caused by the 
work of collective forces of lES. These forces are depen- 
dent on the disks velocities and their work transforms the 
motion energy of lES into the internal energy [6] . 

For generalization of results of researches of hard disks 
it was necessary to pass to the model of nonequilibrium 
systems of potentially interacting of elements. Study- 
ing of disks systems has shown that for this aim it is 
necessary to represent the model of nonequilibrium sys- 
tem as a set of interacting lES. That model possesses 
the big generality. Indeed, the base of statistical physics 
was created by Gibbs using that model [7] . In particular 
it has helped to introduce the distribution function into 
classical mechanics basing on concept of probability of 
beginning conditions for lES. 

The method of Gibbs's ensembles is applicable for 
studying equilibrium systems in the absence of an en- 
ergy exchange between lES. But the exchange of energy 
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between lES in case of nonequilibrium systems is re- 
sponsible for an equilibration. Nevertheless the splitting 
of nonequilibrium system of potentially interacting ele- 
ments into lES allows reducing the proof of irreversibility 
to the proof of existence of non-potential part of collec- 
tive force between lES. It will be proved with the help of 
the equation of interaction of systems (UVS) which was 
obtained directly from the energy conservation law [8] . 

Basing on UVS and the model of nonequilibrium sys- 
tem as a set of lES the new approach to the analysis of 
nonequilibrium systems has been offered. The following 
assumptions and conditions are put into the basis of this 
approach: 1). lES energy should be submitted by the 
sum of internal energy and energy of lES motion as a 
whole; 2). Each element of system should be fixed for 
corresponding lES without dependence from its mixing; 
3). During all process of equilibration all of subsystems 
are considered as equilibrium. 

The first condition is necessary for introduction inter- 
nal energy into the description of dynamics of systems, 
as the parameter which is necessary to correct describ- 
ing an energy exchange between lES. The second con- 
dition allows avoiding a problem of redefinition lES due 
to mixing particles. Last condition is known from ther- 
modynamics. It removes the problems connected to the 
description of system under condition of infringement of 
lES equilibrium. 

Below the substantiation of this approach is submitted. 
The UVS equation and expression for dissipative force 
determining change of internal energy lES are obtained. 
It is shown how basing on UVS and D'Alambert princi- 
ple, the obtaining of the generalized Liouville equations 
for lES is possible. Formulas for entropy of nonequilib- 
rium systems are submitted. Connection between UVS 
and the basic equation of thermodynamics is analyzed. 
Here we also show how and why the offered approach 
has allowed to substantiate thermodynamics within the 
framework of laws of classical mechanics. 

The problem of the description of nonequilibrium sys- 
tems is a key problem for physics as a whole. Therefore 
alongside with mathematical calculations the discussion 
of the reasons why this problem could not be solved on 
the basis of canonical Hamilton equations is offered. 
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II. A SUBSTANTIATION OF THE APPROACH 
IDEA 

The classical mechanics is based on the abstract con- 
cept of the elementary particle and the law of conserva- 
tion of energy. Use of these concepts allows constructing 
classical mechanics on the basis of postulate: "Work of 
forces of reaction is always equal to zero on any virtual 
displacement which is not breaking set kinematics con- 
nections" [9]. Taking into account a condition of con- 
servatism of the active forces leads to a principle of the 
least action, equations of Lagrange, Hamilton and Liou- 
ville [9-12]. These equations describe systems' dynamics 
in equilibrium states and near by. But all attempts to de- 
scribe the dynamics of noncquilibrium systems with their 
help collide with a serious complexities. The analysis of 
these attempts allows assuming that some restrictions 
used at creation of a mathematical formalism of classical 
mechanics are unacceptable for the description of evolu- 
tion of nonequilibrium systems. Our task is to determine 
these restrictions and to discover a way of construction 
of the theory allowing their elimination. 

The evolutionary processes in nonequilibrium system 
are caused by internal forces and streams of energy cre- 
ated by them. These streams is impossible to describe 
with the help of a Hamilton formalism. One of the 
reasons is impossibility of application of a formalism of 
Hamilton for parts of system [9]. There is also other 
not less important reason. The Hamilton formalism is 
constructed on the basis of differential D'Alambert prin- 
ciple under condition of conservatism of collective forces 
[9]. This requirement excludes irreversibility automati- 
cally. But performance of a condition of conservatism 
of forces is strictly proved only for the equilibrium sys- 
tems and in the approach of the theory of disturbances. 
In all other cases including a case of nonequilibrium sys- 
tems, the strict proof of conservatism of collective forces 
is absent [9-12]. Therefore it is necessary to search for 
such approach to the analysis of nonequilibrium systems 
which allows to exclude the using of the requirement of 
potentiality of collective forces. 

It was already mentioned that when the model of sys- 
tem of disks is presented by plurality lES, the irreversibil- 
ity can be shown by analytical way. It is connected 
with non-potential collective forces between lES changing 
their internal energy [5, 6]. It means that for the analy- 
sis of nonequilibrium, systems instead of model of system 
consisting of elementary particles it is necessary to use 
the model of system in the form of plurality of lES. For 
such model the solution of a problem of irreversibility is 
reduced to the proof of presence of non-potential forces 
between lES. 

Necessity of splitting of system on a set of equilibrium 
subsystems at the description of its evolution follows from 
a statistical physics and a physical kinetics. In nonequi- 
librium systems for each physical point there corresponds 
the local velocity of a stream determining energy of a mo- 
tion in this point. Besides in this point the internal en- 



ergy is exist. This energy is determined by energy of the 
chaotic motions of particles. The nature of these two en- 
ergies is various. The first is caused by transport of mass 
in the external field. The second energy is determined by 
a chaotic motion of particles. Therefore the distribution 
function at local equilibrium is determined by energy of 
a motion of lES center of masses (CM) and energy of a 
chaotic motion of particles [13]. The relative motion of 
lES disappears at an establishment of equilibrium. I.e. 
energy of motion of lES is the parameter describing a 
rate of a nonequilibrium. 

It is possible to understand the necessity of splitting 
of lES energy on two types: internal energy and energy 
of CM motion from the analysis of two body system. As 
it is well known the problem of two bodies could not 
be solved directly in a laboratory coordinate system. It 
is related to the nonlinearity arising due to interaction 
between bodies. But the problem can be easily solved by 
transition into CM in which variables are parted. Such 
transition is equivalent to representation of energy in the 
two forms: energies of a motion of system as a whole 
and its internal energy. In the absence of exterior forces 
the motion energy of system is constant and a task is 
reduced to definition of the relative motion of particles. 
The problem can be solved by differentiation on time of 
the energy presented in the form of motion energy and 
an internal energy. 

So, we will build the approach to analyzing of the dy- 
namics of nonequilibrium systems within the frame of 
classical mechanics by the next way. A requirement of 
absolutely elasticity of the elements we will exclude by 
replacement of model of system of elementary particles 
on model of the system consisting from equilibrium struc- 
tured particles, i.e. IBS. Then, having presented of the 
lES energy as the sum of energies of its motion and an 
internal energy, we shall find the UVS. Using UVS the 
equations of Lagrange, Hamilton and Liouville will be ob- 
tained without use of the requirement of conservatism of 
collective forces. 

III. THE SYSTEM MOTION EQUATION 

Let us obtain the motion equation for particle in an 
external field. Let us presume that E = mu^/2 -I- U{r) — 
const is a particle energy. Here m is a mass; T = m.v^/2 
is a kinetic energy; U{r) is a potential energy. Then from 
eq. E = 0, we will obtain: 

v{mv + dU/dr) = (1) 

The eq. (1) is a balance equation of the kinetic and 
potential energies. It is carried out if the condition takes 
place: 

mv = -dU/dr (2) 

It is Newton equation (NE). This equation is integrable 
because a variables were separated. The separation of 
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variables became possible in due to splitting of energy 
into kinetic and potential components, each of which de- 
pends on the different variables. We see that the sum 
of active and inertial forces is equal to zero. The par- 
ticle moves lengthways of a gradient of potential. The 
work on the closed contour is equal to zero. Therefore 
the dynamics of a particle is reversible. 

Let us consider a system consists of A'' potentially in- 
teracting elements; the mass of each clement is equal 
to 1. The force between any two elements is a central 
and determined by the distance between them. Energy 
of system consist of the sum of kinetic energy of ele- 

N 

ments, Tjv = J2 ''nVi^/2, their potential energy in a field 

i=l 

of external forces, - J/jv^"", and the potential energy 

N-l N 

of their interaction UN{rij) = J2 J2 Uij{rij), where 

i=i j=i+i 

= '"i ^ '"i - is a distance between elements i and j. So, 
E = En + J7^"" =Tn + Un + U""^" = const. The time 
derivative of the energy will be as follows: 

TV 

Y.'^iFi = (3) 

i=l 

Where = mvi + E Pij + F?"'' is effective force for i 

N 

particle; f/'^™ = J2 ViF^"''; Fij = dUN/drij; F["^ri) = 

i=l 

The eq. (3) can be treated as orthogonality of the vec- 
tor of effective forces with respect to the vector of veloci- 
ties of elements of the system. If there are no restrictions 
imposed on the Vi directions, the requirement = is 
satisfied [9]. Then from eq. (3) we obtain: 

N 

mvi = -Y,v,Fi-Fr (4) 

i=l 

It is NE for the system's elements in non-homogeneous 
space. It leads to conclusion that the motion of an ele- 
ment of system is determined by the force which equal to 
the sum of vectors of forces, acting from all other parti- 
cles and external force [11, 12]. 

Let us consider the motion of a system as the whole 
in a field of external forces. As well as for an elementary 
particle its kinetic energy is determined by the motion 
of CM. But except of this energy the system possess of 
internal energy. This energy has other nature because it 
is connected with interactions between particles but not 
with the external field of forces. The internal energy as 
well as energy of CM motion should vary due to the work 
of the external forces. Really, the external field will make 
work on change of internal energy due to motion of the 
particles relative to the CM. But velocity of CM motion 
does not depend on motion of elements inside of system. 
In connection with it, the energy of system should be 



presented as the s\im of energy of motion of elements 
relative to the CM and the energy of CM motion. Ac- 
cording with such dividing the system's energy into two 
types, the external forces also will be represented by the 
sum of two forces: the force changing the velocity of CM 
and the force changing the internal energy. 

Let us take into account the equality: Tjv = 

N N-l N 

E mViV2 = MmV^ + {m/N) E E 4 (a), where 

i=l i=l j=i-\-l 

Vjv = -Rjv = 1/-^ E '^j "are velocities of the CM; Rn - 

i=l 

are coordinates of the CM; Vij = fij. 

Let's designate: En = T*[ + E"™" where = 
T^'" + [/at is entrance energy, is a CM kinetic en- 
ergy. The velocity of elements is Vi — Vi + Vn where Vi 
is a velocity of particle relative to the CM. Then: Tjv = 

N N N 

MnV^+tuVn E E 'mvf/2. Because Y,Vi=0 then 

i=l i=l i=l 

N N-l N 

from (a) we have: Y,mv1/2 = 1/{2N) E E ^'y - 

z=l z=l j=i-\-l 

Thus the total kinetic energy of relative motion of par- 
ticles is equal to the sum kinetic energies of their mo- 
tions relative to CM. Because = fij = fi — fj, 
where fi,fj - are coordinates of the elements with re- 
spect to the system's CM then UN{rij) = UN{fij) = 

JV-l N N 

UN{fi) and E E '"ijFijinj) = Y.ViFi{fi), where 

i=l j=i+l i=l 
N 

Fi = dUn/df^ = E dUN/dnj. Thus we have [14]: 

j:^i,j = l 

VnMnVn + E^' = -VjvF^"" - (5) 
Here F^"" = E Fr^R,fi), E^' = E Vi{mii + F{f)i), 

i=l i=l 

N 
i=l 

The eq. (5) represents balance of energy of system in 
a field of external forces. The first term in the left hand 
side determines change of kinetic energy of system. The 
second term determines the change of internal energy of 
system. The first term in the right hand side determines 
the work of forces changing energy of motion of system 
as a whole. The second term determines the work of 
forces changing internal energy of system. This work is 
connected with the motion particles of system relative 
to CM in the external field. Thus the work of external 
forces changes Tj^ and iS]^*. 

When the external forces are absent the eq. (5) will 
split on two independent equations: the first is the equa- 
tion of CM motion and the second is the equation of 
motion of particles relative to the CM. 

Let us note that the eq. (5) can be obtained directly 
basing on the NE for elements. For this purpose we shall 
multiply the eq. (4) on the corresponding velocity. After 
summation the obtained equations for all particles we 
shall have the eq. (5) (if we have summarized the eq. (4) 
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without multiplying it on velocity in this case the internal 
forces in the second term of the eq. (5) will be lost [15]). 
It is confirms validity of the equation (5). 

Let us compare dynamics of a particle and dynamics of 
system. As it follows from the eq. (1, 4) the force acting 
on the particle is potential. The particles dynamics is 
determined by kinetic and potential energies. Obviously 
it is impossible to find particles ' trajectory if not to divide 
the energy into the potential and kinetic parts. 

In agreement with the eq.(5) the work of external forces 
determining motion of system goes both on changes of its 
CM velocity and on change of internal energy. This force 
can be divided on two forces. The first force is potential. 
The momentum of lES is change by this force. The sec- 
ond force changes the internal energy. It is non-potential 
force. Thus it is impossible to describe dynamics of sys- 
tem if not to divide the energy of system on three types: 
kinetic energy of motion of CM, internal energy and sys- 
tem 's potential energy in external field. 

The equation of motion of system in an external field 
can be obtained from the eq. (5). It may be expressed 
as: 

MnVn = -i^""^ - [fi^jj'' + ^'"""Wn/V^ (6) 

Let us call cq. (6) as generalized NE (GNE) for the 
structured particle. The GNE is reduced to the NE if 
one neglects the relative motion of elements, i.e. when 
the internal energy does not change. In this case the 
dynamics of system is similar to reversible dynamics of 
an elementary particle. 

Thus, it is necessary to present the energy of a system 
as a .sum of energy of motion and internal energy for 
determination a system dynamics in an external field. 
The external forces also need to be splitting into the forces 
changing these types of energies accordingly. Otherwise 
variables are not parted. 

IV. THE EQUATION OF INTERACTION OF 
TWO SUBSYSTEMS 

As it follows from the kinetic equations the dynamics 

properties of the nonequilibrium systems' arc determined 
by the energy of motion and internal energy in each phys- 
ical points of the system [7, 13]. These parameters can be 
entered into the description of dynamics of system if the 
system is presented by a set of lES or, in other words, by 
a set of the structural particles. Using model of system as 
a set of lES help us "select" nonequilibrium effects into 
interactions of lES. These interactions can be determined 
by the eq. (5) if external forces will be replaced by the 
forces between lES. Let us show how it can be done. 

Let suppose the system consists of two lES - L and K. 
We take all elements to be identical and have the same 
weight 1, and L to be a number of elements in L - lES, 
K -is a number of elements in K -lES, i.e. L -\- K = N, 

L K 

Vl = l/L^Vi and Vk = K ^ Vi - are lES's velocities 

i=l i=l 



with respect to the CM of system. The velocity of the 
system's CM we take equal to zero, i.e. LVl + KVk = 0. 

We can represent the energy of the system as = 
El + Ek + t^*"* = canst, where El and Ek are the 
lES energy, and C/™* - is the energy of their interac- 
tion. According to the eq. (6), the energy of each 
lES can be represented as El = T^l!' + ^Y^' ^ = 
Tf^ + E'^', where Tf - MlVI/2, Tf^ = 
Ml = rnL, Mk = niK. iJ*""- is the internal energy 
of a lES. The i?™* consists of the kinetic energy of mo- 
tion of the elements with respect to the CM of lES - 
T'"^ and their potential energy - C/*"^ i.e. = 

Tins + U'ns^ ^^^^^ = i: ^ Ui,j,{n,j,), 

Ur = e\ E U.^jKiri^jK)- The energy is 

K L 

determined as f/'"* = E E ^jLjKifjLjK)- Indexes 

jK = l jL=l 

.jk, .jLjiKTiL determine belonging of the elements to cor- 
responding lES. In equilibrium we have: T*'' = 0. Hence, 
if the system aspirates to equilibrium, then T*'' energy for 
each lES will be transformed into the internal energy of 
lES. 

We can obtain the equations of dynamics of L and K 
of lES in the following way. Let us differentiate energy 
of system on time. In order to find the equation for L - 
lES, at the left hand side of obtained equality we have 
kept only those terms which determine the change of the 
kinetic and potential energies of interaction of elements of 
L - lES. We replaced all other terms in the right hand side 
and combined the groups of terms in such a way when 
each group contains the terms with identical velocities. 
In accordance with NE (sec eq. (5), the groups which 
contain terms with velocities of the elements from K- 
lES are equal to zero. As a result the right hand side of 
the equation will contain only the terms which determine 
the interaction of the elements L-IES with the elements 
ii'-IES. The equation for if-IES can be obtained in the 
same way. Then we execute replacement of variables: 

= Vi -\- V and take into account equality (a). As a 
result we will have [8, 14]: 

VlMlVl + Er = -^L - Vl* (7) 

VkMkVk + Er = ^K + (8) 
Here * = E $l = E ^^i.F^; = 

K K L 

E ^iK-^lK^ ^IL ~ E ^iLjK'l ~ E PiLjKI 

iK=l jK=l ii, = l 

L-1 L 

■pins 1,. . [""iiLiE. I 

^L ~ ^ilJl[ Z I" 

+ F..,. ] ; Er = e' E vi^,^ [^%^ + 

iK = '^ jK=iK+l 

+ FiKjK]- 
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The eqs. (7, 8) arc the UVS. They describe energy ex- 
change between lES. Independent variables of UVS are 
macroparameters, coordinates and velocities of lES mo- 
tion, and also microparameters, coordinates and veloci- 
ties of elements. UVS binds together two types of the 
description: at a macrolevel and at a microlevel. The 
description at a macrolevel determines of lES dynamics 
as a whole and at a microlevel determines dynamics of 
lES elements . 

The potential force, "if, determines motion of lES as a 
whole. This force is the sum of the potential forces acting 
on elements of one lES at the side of another lES. 

The forces which determined by the terms $l and 
^K, transformed of the motion energy of lES to inter- 
nal energy as a result of chaotic motion of elements of 
one lES in a field of the forces of another lES.They are 
non-potential forces which dependent on velocities and 
cannot be expressed as a gradient from any scalar func- 
tion. These forces are equivalents to dissipative forces. 

The motion eqs. for lES can be obtain from eqs. (7,8): 

MlVl = - ulVl (9) 

MkVk = * + axVK (10) 

where = (S™^ + <i>L)/Vl, uk = {^k - E'D/Vk^ 

The eqs. (9, 10) are GNE for interacting of systems. 
The second terms in the right hand side of the equa- 
tions determine the forces changing internal energy of 
lES. These forces are equivalent to friction force. Its 
work consists of works of chaotic motion of elements one 
lES in a field of forces of another lES. The efficiency of 
transformation of energy of relative motion of lES into in- 
ternal energy are determined by the factors "ai,", "qk"- 
If the relative velocities of lES elements are equal to zero 
the force of friction is also equal to zero. 

For the description of dynamics of nonequilibrium sys- 
tems we instead of traditional model of system have taken 
model of system as a set of interacting lES. Here a role 
of particles is carried out by lES. Their motion is deter- 
mined by the eqs. (9, 10). 

The state of system from a set of lES can be deter- 
mined by the point in the phase space which consists 
of 6i? — 1 coordinates and momentum of lES, where R 
is a number of lES. Let us call this space as S'-space 
to distinguish it from usual phase space for elementary 
particles. Unlike usual phase space the S'-space is com- 
pressible though total energy of all elements is a constant. 
The rate of compression of 5-space is determined by the 
rate of transformation of motion energy of the lES into 
their internal energy. Thus the volume of compression of 
5-space is determined by energy of the lES motion. 

Let us compare GNE and NE. The NE follows from 
representation of energy of an element as the sum of po- 
tential and kinetic energies. The transformation of one 
energy type to another is determined only by the poten- 
tial force under condition of conservation of full energy. 



The GNE determines motion of lES not only at trans- 
formation of kinetic and potential energies of lES, but 
also at transformation of lES internal energy. These 
transformations are determined by two forces. The po- 
tential part of force changes velocity of lES's CM. The 
non-potential part of force changes its internal energy. 
Thus, for the structured particle external force will con- 
sist of two parts: potential and non-potential. Each of 
them, has collective character. 

The GNE passes into NE when the relative velocities 
of elements lES can be neglected. It is also have a place 
when distances between lES are big enough [16]. 



V. THE GENERALS OF LAGRANGE, 
HAMILTON AND LIOUVILLE EQUATIONS FOR 

lES 

The dynamics of the nonequilibrium system submit- 
ted by a set of lES will be determined by generalized 

Lagrange, Hamilton and Liouville equations for lES. Be- 
fore explaining a way of obtain of these equations, we 
shall remind under what conditions the canonical prin- 
ciple of Hamilton for system of elementary particles was 
deduced [9]. 

According to the differential principle of D'Alambert 
"the work of effective forces including inertial and active 
forces, at virtual motions elements of system is equal to 
zero for all reversible virtual motions compatible with re- 
strictions on dynamics" [9]. From D'Alambert principle 
with the help of NE the integrated principle of Hamilton 
is deduced. For this purpose the integral on time of the 
virtual work made in system by effective forces is equated 
to zero. Integration on time is carried out provided that 
external forces possess power function. It means that the 
canonical principle of Hamilton is fair only for cases when 
^ FiSRi = —6U where i is a number of particles, and Fj 
- is a force acting on this particle [9]. But for lES we 
have no right to demand a carrying-out of a condition of 
conservatism of its forces. From the mathematical point 
of view it is because the Fi^Ri cannot be equal a full 
differential for lES. Prom the physical point of view it 
is due to the non-potential force which responds to the 
changes of the internal energy of lES. Therefore obtaining 
of the Hamilton equations should be carried out basing 
on GNE. 

The principle D'Alambert for lES sounds so: the sum 
of works of all forces of interaction IBS at their virtual 
motions compatible with restrictions on dynamics and in 
view of change of internal energy of IBS is equal to zero. 

In according to the eqs. (9, 10) and Liouville formula 
from [6, 9], the generalized Liouville equations in 5-space 
can be written as: 

df/dt = -fdF/dV (11) 

Here /- is a distribution function for lES, F -is a non- 
potential part of force acting on lES. 
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The right hand side of the eq. (11) is not equal to 
zero as forces of lES interaction which transformed of 
the energy of its relative motion in its internal energy, 
depend on velocities of elements. The S'-space which de- 
termined by the coordinates and velocities of lES will be 
compressed. Thus system of lES will converge to equi- 
librium. 

The impossibility of return of internal energy of lES in 
its energy of motion is caused by impossibility of change 
of lES motion due to motion of its elements. At an es- 
tablishment of equilibrium the motion energy of lES will 
disappear. In this case the description in 5-spacc coin- 
cides with the traditional description in usual phase space 
for elements of system. 

VI. UVS AND THERMODYNAMICS 

It is possible to come to the idea about necessity of re- 
placement of model of elementary particles by the model 
of interacting lES from the analysis of the basic equation 
of thermodynamics. According to this equation the work 
of external forces acting on the system are splitting on 
two parts. The first part is connected to reversible work. 
It can be associated with change of energy of motion of 
system as a whole. Second part of energy will go on 
heating. It is connected to internal degrees of freedom of 
system. To this part of energy there corresponds internal 
energy of lES. Here we would like to show, how basing 
on the UVS it is possible to come to the thermodynamics 
from the classical mechanics. 

Let us take the motionless nonequilibrium system con- 
sisting of "i?" of lES. Each of lES consists of enough 
plenty of elements. Let dE is a work which has been 
done at system. In thermodynamics that term is an in- 
ternal energy of a system (do not confuse E with the iJ'"* 
- internal energy of lES). The dE is determined by the 
basic equation of thermodynamics as: dE = dQ — PdY 
[8]. Here, according to common terminology, E is energy 
of a system; Q is thermal energy; P is pressure; Y is 
volume. 

As well as the basic equation of thermodynamics, UVS 
also is differential of two types of energy. According to 
the UVS the volume dE is redistributed inside of the 
system in such a way that that one its part goes on change 
of energy of relative motion of lES, and another part 
changes its internal energy. The first term in the left 
hand side of UVS is a change of kinetic energy of motion 
of a lES as a whole, dT*^. This term corresponds to the 
term PdY. Really , dT*'' = VdV = VVdt = Vdr = PdY 
[7]. 

If the potential energy is a homogeneous func- 
tion of a second order of the radius-vectors, then 
in agreement with the virial theorem [10], we have: 
^ins ^ 2f^"' = 2W'^\ The line denotes the time aver- 
age. 

Let us consider the system near to equilibrium 
state. The average energy of each element is 



^tns ^ E'^'^'/N ^ kT^"'' where A^ is a number of el- 
ements. As the increasing of the internal energy is 
determiucxl by the volume dQ, then we will have: 
«To""[d.E'^"Vr^"^] ~ To"«[d^;o/^^o], where vq is the 
average velocity of an element, and dvo is its change. For 
the system in the closed volume we have: dvo/vo ~ dT/T, 
where F is the phase volume of a system, dV will 
increase due to increasing of the system's energy on 
the value, dQ. By keeping the terms of the first or- 
der we get: dQ « T^'^'dT/T = Tg^^dlnF. By definition 
dlnV = where 5'"* is an entropy [8]. So, near to 

equilibrium state we have d£'"'* = dQ ^ r(5"*d5"". 

The entropy increasing, AS*, for nonequilibrium system 
is completely determined by the energy T**" passing into 
Therefore AS' can be determined by the formula 

[7, 14]: 

^s=j:{mf:lY.^d^} (12) 

L=l /s=l s 

Here E^ is the kinetic energy of L-IES; Nl is the num- 
ber of elements in L-IES; L = 1, 2, 2>...R; R is the number 
of lES; s is a number of the external element which in- 
teraction with clement k belonging to the L-IES; L"^ is 
a force, acted on fc-element; Vk -is a velocity of the k- 
element. 

In agreement with eq. (12), the entropy is determined 
by the energy of the relative motion of lES, transformed 
in an internal energy as a result of work of non-potential 
part of forces between lES. 

To obtain equation for the entropy production we 
take into account that: AS = AQ/T. Thus dS/dt = 
[dE^"'^ /dt]/{kE''"'^), where fc is a coeflScient. It is possi- 
ble to express this formula through the work of forces of 
interaction lES. Let Eq is a full system's energy, E^^^ 
is a beginning energy of relative motion of L-IES. In 
accordance with eqs. (10,11) the rate of increasing of 
the internal system's energy is equal to: ( = ^i^—i^L- 
The internal energy of a system is equal to: LJ'"* = 
Eq — where a sum of lES en- 

ergy of relative motion. But Y^^=i = Co — Jo C{t)dt, 
where Co = J2l=i ^lo- Then entropy production for the 
system, Qprod = dS/dt, can be write: 

Qprod = D/{1 -Do + JoD{t)dt), where D = C/Eq, 
Do = Co/E. 

The energy of T**" characterizes the rate of non- 
equilibrium system while LJ'"* characterizes its degree 
of equilibrium. 

For maintenance of nonequilibrium system in a sta- 
tionary state, loss of energy L^*'' must be compensated 
by the inflow of external energy. It is possible to cre- 
ate with the help of contact of system with heater and 
due to outflow of radiation. Then the stationary state of 
nonequilibrium system is characterized by the formula: 
Qprod = l^"-! — I £•-1-1) where Q- is entropy outflow, q+ -is 
inflow of entropy, Qprod is entropy production, determined 
by the formula (12). 
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The process of transformation E^"" to the energy of 
lES motion is forbidden by the law of conservation of 
lES momentum [10]. Really, momentum of lES cannot be 
changed by the forces which change the internal energy. 
Therefore the work along the closed contour for lES not 
equal to zero. As a result the concept of entropy as the 
rate of increasing of lES internal energy due to energy 
of its relative motion will appear within the framework of 
classical mechanics. It is impossible to come to entropy 
concept in the frame of classical mechanics if base on the 
model of system as a set of elementary particles because 
for such model the determining of E^"^ is impossible. 

VII. CONCLUSION 

To connect thermodynamics with classical mechanics 
it is necessary to refuse from the idealization of an ele- 
mentary particle and a condition of conservatism of col- 
lective forces. It can be made by replacement of system 
of elementary particles on a system of lES. 

Collective forces between lES are consisting of two 
parts. The potential part of force changes energy of lES 
motion. The non-potential part of force changes inter- 
nal energy. The change of internal energy is caused by 
chaotic motion of particles of one lES in a field of forces 
of others lES. Therefore at equality to zero of velocities 
of particles in relative to CM of lES, the changes of its 
internal energy do not occur. 

The Lagrange, Hamilton and Liouville equations on 



the basis of UVS are deduced. These equations describe 
dynamics of the nonequilibrium system in the compress- 
ible S'-space where points are corresponds to the mo- 
mentum and coordinates of the lES. We call this space 
as 5-space in accordance with the entropy designation. 
Compressing of ^-space is explained by the transition of 
energy of lES motion into their internal energy. 

The following mechanism of irreversibility can be pro- 
posed: the energy of relative motion of lES is transformed 
into their internal energy as a result of the work of the 
non-potential part of a force between of lES. The system 
equilibrates when relative motion of lES disappears. 

The explanation of the First law of thermodynamics is 
based on the fact that the work of subsystems' interaction 
forces changes both the energy of their motion and their 
internal energy. 

The explanation of the Second law of thermodynamics 
is based on the condition of irreversible transformation of 
the subsystems' relative motion energy into their internal 
energy. 

Thus within the limits of the accepted restrictions the 
offered approach has allowed to expand a formalism of 
classical mechanics on nonequilibrium systems and to 
connect it with thermodynamics. 

Difficulties on a way of the solution of a problem of 
irreversibility are resolved at refusal of the assumption of 
absolute elasticity of particles and potentiality of collec- 
tive forces between lES. 
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